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FORECASTER, A MARKOVIAN MODEL TO ANALYZE 
THE DISTRIBUTION OF NAVAL OFFICERS 

by 

Paul R. Milch 

1. INTRODUCTION AND BACKGROUND 

In an earlier report [4] the author constructed an analytical model based 
on Semi-Markov Processes representing the career paths of military officers 
in general and Surface Warfare Officers (SWO's) of the U.S. Navy in 
particular. The model was based on a rectangular grid representation of the 
entire career structure of a group of officers. In the rectangular grid, rows 
represent the various activity types (e.g., sea duties and shore duties) in 
which an officer of the group whose career structure is modeled may be 
engaged from time to time. Columns of the grid stand for the succession of 
tours that makes up the career of any military officer. In this way, an officer's 
career may be represented by a path through a number of nodes (billets) 
whose first coordinates show the activity types and the second coordinates the 
tour numbers. 

Associated with each node on the grid is a specific non-negative integer, 
called the tourlength, which determines the length of time the officer spends 
in that node or billet. A zero tourlength means that the node or billet is 
infeasible. This way, a matrix of size A by R determines the tourlength matrix 
where A and R are, respectively, the number of activity types and the number 
of tours in the entire career structure. A similar A by R matrix may represent 
the number of officers, called incumbents, currently occupying billets in all 
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the feasible nodes. Yet a third A by R matrix of non-negative integers 
provides the numbers of accessions or recruits entering the system at future 
times. Typically, this latter matrix is all zeros, except for some elements of the 
first column, because new entries into the system should start with the first 
tour. However, in some officer commimities, such as the Navy Nurse corps, 
for example, there are significant numbers of "recalls" and other "lateral 
entries" who join the community in later tours. Finally, the progression of 
advancement on the career grid is governed by a sequence of transition 
probability matrices which prescribe the probability of an officer currently in a 
specific node (i.e., activity-tour combination) transiting to another (possibly 
the same) activity type upon finishing his/her current tour. These matrices 
are all A by A and there are one fewer of them than the total number of tours 
(i.e., R-1). A more detailed description of the system is given in Milch [4]. 

In the same report the mathematical details and formulas are also 
described for the expected numbers of officers occupying each node (billet) on 
the grid at some future point in time. This is based on first computing 
probabilities that an officer occupying a node (billet) now will be anywhere on 
the grid t time units later. This computation is based on a Semi-Markov 
Process representation of the position occupied by an officer at time t. When 
the appropriate formulas were programmed for computation in APL on the 
IBM 3033 mainframe computer at the NPS it was found that the computation 
required either excessively large amounts of memory or huge amounts of 
CPU time. This occurred even after various time-saving measures were 
effected. Since it was desirable to make the model available for personal 
computers as well, the necessity of revising the computational core became 
even more acute. 
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Simultaneously with efforts to revise the mathematical core of the model 
two attempts have been made to evaluate the practical usefulness of the 
model apart from its computational speed (or rather lack of it). Thus Johnson 
[3] built an interactive user friendly interface for the model now labeled 
FORECASTER, thereby greatly enhancing its usability. He also demonstrated 
how FORECASTER may be used by analyzing the currently topical issue of 
joint duty assignments effected by the recently passed Goldwater-Nichols 
Department of Defense Reorganization Act. Title IV of this act mandates the 
creation of so-called joint duty officers in all four military services and 
prescribes their professional education, as well as their use by the services. 
Johnson examined the effect of the new law on the SWO community: its 
ability to fill its "fair share" of joint duty billets and how that may affect the 
traditional SWO career path. More recently, Drescher [1] examined the same 
problem in another Navy officer community, namely that of Tactical 
Aviation pilots and flight officers. At the same time, Drescher reorganized 
the APL functions used in FORECASTER and provided a more thorough 
documentation for the model. 

Concurrently with Drescher's work, this author has revised the 
mathematical core of the model. This new mathematical formulation, a 
conceptually much simpler one than the original presented in Milch [4] 
turned out to be significantly faster computationally as well and that is the 
subject of this report. It is also hoped that the new mathematical formulation 
of the model will make it easier to effect another significant improvement in 
the usefulness of the model. This issue has to do with the way the user must 
input data into the current model, as well as the way the user must make use 
of the output results. 
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Briefly summarized, this aspect is due to the fact that the model is based 
on a sequence of tours that the officers undergo. This is a realistic simulation 
of what actually happens, but it also means that incumbents data (i.e., 
numbers of officers currently occupying billets) must be categorized by 
activities and tours. Unfortunately, the usual method of recordkeeping of 
officers in any Navy community is not organized by tours, but rather by grade 
and years of commissioned service. Likewise, results of the model categorize 
officers distributed at some future time by activity and tour number. These 
results must then be compared to billet requirements that are usually not 
organized by tour number but by grade. 

Johnson [3] and Drescher [1] demonstrated that the data is available by 
tour number from the source where all officer data originates, the Officer 
Master File. However, the fact remains that standard personnel planning by 
all Navy offices where such work is done is based on officer and billet data 
that are not organized by tour numbers and this practice is not likely to 
change in the near future. For this reason it is desirable to convert the model 
to accept input and produce output by grade and/or years of service rather 
than tour number in order to accommodate user needs. 

This effort is currently under way. In the meantime, the new 
mathematical core was designed to accept data as the original model did 
which is described in Milch [4] and outlined at the beginning of this section. 
The next section describes the mathematical details of the model according to 
the newly created APL code available on either the IBM 3033 mainframe 
computer or a personal computer. 
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2. DATA PREPARATION 

As a point of departure, the notation of the author's previous report [4] 
will be adopted wherever appropriate. Thus, the tourlength matrix, L, has the 
elements 

_ duration of time spent in activity type i 
^in - by an officer on his/her n‘*' tour, 

for i = 1, ..., A and n = 1, ..., R. 

The incumbents matrix, Q, is defined as having elements 

number of officers at time zero occupying 
*1'"“ activity type i billets during their n‘^ tour, 

for i = 1, ..., A and n = 1, R. 

Further, the accessions matrix, C, is defined as consisting of elements 

number of officers accessed during a time period 
“ directly into an activity type i and tour number n, 

for i = 1, ..., A and n = 1, ..., R. 

Most of the time, matrix C will have positive elements only in its first 
column (for n = 1), since a newly accessed officer must usually start in tour 
number one. It is interesting to note, however, that this is not always the 
case. For example, officers are sometimes recalled shortly after their 
separation from the service or are laterally transferred from one community 
to another and may thus join the community being analyzed in a tour other 
than their first one. 

Finally, the transition probability matrices, P(n) are defined to have 
elements. 
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probability that an officer transfers to an activity type j 
Pij(n) = on his/her n+l®‘ tour upon completion of his/her 
n‘*^ tour in an activity type i, 

for i,j = 1 , A and n = 1, R-1. 

The construction of the model will be described in steps like those of an 
algorithm as that is most suitable for the procedure which turned out to be 
both simple and computationally efficient. 

Step 1: Unravel the tourlength matrix, L, by columns, starting with the 

A 

first column, followed by the second, etc., into a vector £ . More 
precisely, define 

A A A 

£ = (^1, ...,^k) 

where 

A 

/ k = ^in when k = (n-l)A+i 

for i = 1, ..., A and n = 1, ..., R-1 so that k = 1, 2, ..., K = AR 

Step 2: Define the indices ki, k 2 , ..., kK» as the first, second, ..., last values 

A 

of k for which ^ k is positive. These are called the "ascending 
ladder indices or epochs" by Feller [2]. Then the abbreviated 

A 

version of £ may be defined by omitting its zero components. In 
other words, 

where 



R A 

K* = X X Sin with 
n=l i=l 




if /in>0 

if 



It will be useful to define the partial sums of the components of 
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t* = 





as well. Let 

U U 

So = 0 and Su = X =X ^ L for u = 1, ..., K* 

^ A j=i I j=i ^ 

since t. = I k- Feller's terminology the Su are the "ascending ladder 
heights." 

Step 3: Next the vector is expanded by repeating each of its 

components exactly as many times as indicated by its value. That 
is, if = (3, 1, 5) then the expanded vector will have nine 
components with the first three equal to 3, the fourth equal to 1 
and the last five equal to 5. In general. 



t* = {I 



1 




) 



where 



K* = X X An 

n=l i=l 

The components of may formally be defined as 



C* = C when Su-i + 1 < j< Su and u = 1, ..., K*. 

j K 



For example, when u = 1, ^ = ( for all j = 1, ..., Si = ^ . This means that 

j '^1 



for the above example i =3, for j = 1, 2, 3. 

j 

Step 4: Next a similar procedure is followed with respect to the 

incumbents matrix, Q. First, Q is unraveled by columns into the 

A 

vector q. That is. 
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where 



g=(qi q^) 



A 

qic “ qin ''^hen k = (n-l)A+i 

for i = 1 , A and n = 1, R so that k = 1, K = AR. 

Then g is abbreviated by the omission of those of its components, qj^, for 

A 

which the corresponding / is zero. That is, 

r = (<. •••' 

^ K 

* A 

where q = for j = 1, ..., K*. 

j i A A 

Note that the decision to keep or omit a component, q^, of q is based not 

A A ^ 

on whether qj.>0 or qk=0, but on whether /k>0 or ^k=0- 

Finally, the vector q* is expanded to repeat each of its components as 
many times as indicated by the value of the corresponding component of £*. 
More formally. 



q** = (q 



>W- 
1 ' 



•••/ 



JW- 







) 



where 

q”^ = q* when Su-i + 1 ^j^Su and u = 1, ..., K*. 
j K 

Step 5: Next two "flow rate" vectors are defined by dividing the number 

of incumbents by their respective tourlengths: 

r* = (r*, ..., r* ) and r- = (r"*, ..., r*" ) 

1 K* 1 K** 
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where 



and 



5t* 

r 

j 



* 




t 

j 



j = i. 



K* 



Jt-st- 

q. 

r*’^=— for j = 1, ..., K**. 
j 

j 



Step 6: 



A slightly different procedure is followed with the accessions 
matrix, C. Again, C is first unraveled by columns into the vector 



where 



Ck = Qn when k = (n-l)A+i 

for i = 1, ..., A and n = 1, ..., R so that k = 1, ..., K = AR. 

A A 

Then c is abbreviated by the omission of those of its components, c ^ for 

A 

which the corresponding ^ k is zero. That is. 



c* = (c*, ..., c* ) 
1 K* 



c* = Ckj for j = l, ..., K*. 



) 



where 

Then c”^ is expanded to have K** components, but instead of repeating its 
components (as was the case in expanding into the extra components 
in c** are simply zeroes. More precisely. 



c** = (c , ..., c ) 
1 K** 
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where 




when j = Su-i + 1 and u = 1, 
otherwise. 



K* 



These six steps accomplish the task of rearranging the data originally 
organized for user convenience into formats that fit the need of efficient 
storage and computation. Namely, the tourlength matrix, L, was 
manipulated first to eliminate infeasible billet types and then to repeat its 
positive components in preparation of the computation of the flow rate 
vectors r* and r**. Similarly, the incumbents matrix was first transformed 
into a vector in which any incumbents occupying infeasible billets (due to a 
data recording error or for any other reasons) have been omitted. The 
practical effect of this is that such personnel will be attrited by the model 
during the forecasting process which is probably the appropriate 
representation of what happens to such personnel, anyway. 

The purpose of dividing the number of incumbents by the corresponding 
tourlength and then repeating that quotient as many times as indicated by the 
value of the tourlength is to distribute personnel occupying a billet type 
evenly in terms of longevity (or experience) in that billet. This may not be 
the same as the actual distribution in reality, but the latter may either not be 
available or would present too much trouble for the user to enter into the 
model for the additional accuracy that would be gained thereby. 

The accessions matrix is similarly treated except that instead of repeating 
its components all accessions are placed into the first (i.e., zero) longevity (or 
experience) category and the additional components of c** are simply zeroes. 
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The next task is to consolidate the transition probability matrices, P(n), n = 
1, R-1, into a single large transition probability matrix that will affect the 

transitions of personnel occurring each unit of time. Conceptually the 
simplest way of accomplishing this would be to construct a Markov Chain 
with states corresponding to the components of the vector r*^ Then 
each unit of time components of r** would be transferred in accordance with 
a K** by K** transition probability matrix which would contain all the 
positive elements of all the P(n) matrices, but would otherwise consist of 
mostly zeroes and many ones, representing the fact that personnel cannot 
move backward (or "sideways") in time and once placed in a billet personnel 
must, by necessity, move up by one longevity (experience) level per unit of 
time until their experience level equals the tourlength in that billet at which 
point the appropriate elements of the relevant P(n) matrix would distribute 
them into various activities of their next tour. 

The problem with this approach is that being the sum of all 

tourlengths, is in the range of about 300 or more for most problems for which 
this model is being designed when tourlengths are given in terms of quarters 
(3 months) which is probably the preferred unit of time, with the alternative 
being months (which in turn would further increase the number of states 
three-fold). When forecasting several units of time into the future having 
the number of states in the neighborhood of 300 would tax the computational 
capabilities of most micro-computers in terms of memory and/or CPU time 
and perhaps even of most commonly available mainframe computers for 
which this model could be designed. Therefore, practicality dictates a better 
approach which is less consumptive of CPU time as well as memory. 
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The problem of too high dimensionality will be solved by using the 
Markov Chain approach not to forecast the vector r** but the much lower 
dimensional vector r.*. This reduces the dimensions from K** 
(approximately 300) to K* (approximately 50 or 60). This requires the 
transition probability matrix to be of only K* by K*. 

It is this matrix that is constructed next from the sequence of transition 
probability matrices P(n), n = 1, R-1. This is accomplished in several steps. 

First the matrices, P(n), are made consistent with the information 
contained in the tourlength matrix, L. Namely, if the element hn is zero, 
indicating that activity type i for tour number n is not feasible then the i**' row 
of the matrix P(n) is superfluous and in fact should contain all zeroes. If 
some elements, Pjj(n), are positive for some j = 1, ..., A it must be due to some 
error, most likely an omission to change these elements to zero by a user of 
the model when a formerly feasible billet was made infeasible. At cmy rate, 
^in = 0 implies that Pjjfn) should be zero for all j = 1, ..., A and therefore these 
elements occupy redundant space in P(n). Similarly, the i*^ column of the 
matrix P(n-l) should contain all zero elements and is likewise redundant. 
This results in the following 

A 



Step 1: Define P(n), for n = 1, ..., R-1, to have elements 



Pij(n) = 



p..(n) if 1. i. >0 
in jn+l 



0 



if •^in^jn+l ~ 0 



for i,j = 1, ..., A. 



A 

Then define the matrix P*(n), for n = 1, ..., R-1, to be the same as P(n) except 
that any all-zero row or all-zero column is omitted. Note that the matrices, 
P*(n), may no longer be square matrices and in fact it will be useful to 
introduce notation for their row and column dimensions. Namely, let 
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dr(n) = no. of rows of P’^(n) 
dc(n) = no. of columns of P*(n) 

for n = 1, ..., R-1. 

Also, let dr(0) = dc(0) = 0. 

Further, define the partial sums of these row and column dimensions: 

Dr(n) = £ dr(m) 

m=0 

n 

Dc(n) = ^ dc(m) 

m=0 

for n = 0, 1, ..., R-1. 

Note that 

A A 

dr(n) = X 5in and dc(n) = X Sin+i 
i=l i=l 

for n = 1, ..., R-1. 

Therefore, dc(n) = dr(n+l) or the column dimension of P*(n) is the same 
as the row dimension of P*(n+1). 

Also, note that 

R A R-1 A 

S Sin = X dr(n) + X SiR = Dr(R-l)+dc(R-l) 

n=l i=l n=l i=l 

and similarly, 

R-1 A A R-1 

K* = X X 5in+i =X Sil + X dc(n) = dr(l)+Dc(R-l) 

n=0 i=l i=l n=l 

Therefore, 

Dr(R-l) + dc(R-l) = dr(l) + Dc(R-l) 
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or 

Dr(R-l) - dr(l) = Dc(R-l) - dc(R-l) 

That is, the sum of all but the first row dimensions equals the sum of all but 
the last column dimensions. 

Step 2: The K* by K* dimensional transition probability matrix, P*, may 

now be introduced with all elements zero except the following: 

p* = p*(n) 

^Dr(n-l)+i, dr(l)+Dc(n-l)+j ij 

for i = 1, ..., dr(n) and j = 1, ..., dc(n) and n = l, ..., R-1. 

Note that by necessity the first dr(l) columns of the matrix, P*, must be all 
zeroes. Similarly, the last dc(R-l) rows of the matrix P*, are also all zeroes. It 
is also clear that P* is an upper triangular matrix, in fact all its main diagonal 
elements are also zero. 

The following example may serve as an illustration of the entire 
procedure described above. The three matrices, L, Q and C are given as 





r 


2 


0 


3 


0 


A 






r 


50 




0 


60 


0 




L = 




0 


1 


0 


2 






Q = 




0 


10 


8 


50 






1 


3 


2 


0 


4 


) 








90 


30 


20 


80 


J 










/ 


70 


0 


10 


0 


















C = 




10 


0 


5 


0 




















V 


50 


0 


0 


0 


J 











with i = 3 and R = 4. 

Then the unraveling process produces the K=12 dimensional vectors 

A 

I =( 2 , 0 , 3 , 0 , 1 , 2 , 3 , 0 , 0 , 0 , 2 , 4 ) 
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g = (50, 0, 90, 0,10, 30, 60, 8, 20, 0, 50, 80) 
c = (70,10, 50, 0, 0, 0, 10, 5, 0, 0, 0, 0) 

A 

Examining t shows that the "ascending ladder epochs" are: 

1<i = L1<2 = 3, k3 = 5, k4 = 6, ks = 7, k^ = ll, k7 = 12. 

Therefore, the abbreviated vectors are K* = 7 dimensional, namely: 

r = (2, 3, 1, 2, 3, 2, 4) 

g* = (50, 90, 10, 30, 60, 50, 80) 

c* = (70, 50, 0, 0, 10, 0, 0) 

The "ascending ladder heights" are: 

so = 0. Si =2, S 2 = 5, S 3 = 6, S 4 = 8, S 5 = ll, S6 = 13, S 7 = 17. 

The vectors with repeated components are K** =17 dimensional, namely: 
r* = (2, 2, 3, 3, 3, 1, 2, 2, 3, 3, 3, 2, 2, 4, 4, 4, 4) 
g»» = (50, 50, 90, 90, 90, 10, 30, 30, 60, 60, 60, 50, 50, 80, 80, 80, 80) 

= (70, 0, 50, 0, 0, 0, 0, 0, 10, 0, 0, 0, 0, 0, 0, 0, 0) 

Next, the two flow rate vectors of dimensions K* = 7 and K** = 17 are 
computed as: 

t = (25, 30, 10, 15, 20, 25, 20) 

t* = (25, 25, 30, 30, 30, 10, 15, 15, 20, 20, 20, 25, 25, 20, 20, 20, 20) 

Assume that the following three P(n), n = 1, 2, 3, matrices are also given: 





f 


0 


.5 


.4 


>1 




^ .8 


0 


.1 






7 .2 .6 


.1 




P(l) = 




.1 


.7 


0 




P(2) = 


.8 


0 


0 




P(3) = 


.2 .5 


.2 






1 


.0 


.8 


.1 


J 




1-7 


.1 


0 


J 




1-8 -1 


0 


J 
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The first step is to replace any positive pij(n)'s with zeroes if either ^in = 0 or 



^jn— 0- 





^ 0 


.5 


.4 






^ 0 


0 


0 


>1 




< 0 


.6 


.1 


A 


A 










A 










A 










P(l) = 


0 


0 


0 




P(2) = 


.8 


0 


0 




P(l) = 


0 


0 


0 






lo 


.8 


.1 


J 




1-7 


0 


0 


J 




lo 


0 


0 


J 



Next, all-zero rows and all-zero columns are eliminated: 



P*(l) = 



^ .5 .4 ^ r.8^ 



.8 .1 



P’^(2) = 



.7, 



P’^(3) = (.6 .1) 



Note that 

dr(l) = 2, dr(2) = 2, dr(3) = 1 and dc(l) = 2, dc(2) = 1, dc(3) = 2. 
Therefore, 



Dr(0) = 0, Dr(l) = 2, Dr(2) = 4, Dr(3) = 5 



and 



Dc(0) = 0, Dc(l) = 2, Dc(2) = 3, Dc(3) = 5. 



The K*=7 dimensional square matrix P* is 



P* = 



V 



0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 



/.5 A\ 0 0 

V.8 .1/ 0 0 

0 0 /". 8 \ 0 

0 0 V.7/ 0 

0 0 0 (.6 

0 0 0 0 

0 0 0 0 



0 ^ 
0 
0 
0 



. 1 ) 



0 



Note the way the matrices, P*(l), P*(2) and P*(3), are embedded in the 
large P* matrix, in the form of a "descending staircase" which is typical of the 
structure of the P* matrix. 



3. FORECASTING 

Forecasting will be accomplished in several steps using both K*- and K**- 
dimensional vectors. 
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The first five of these steps are taken for a single value of t, starting with 
t = 0 and then repeated for t = 1, T-1 where T is the forecasting horizon 
for which the distribution of officers is required. After the first five steps 
have been successfully repeated for t = 0, 1, T-1 the actual forecasting 

procedure is complete. Steps 6, 7, and 8 then serve the purpose of rearranging 
the result in the standard A by R matrix format. 

Step 1: First the K*-dimensional vector r*(t) is established, made up of 

those elements of r**(t) that stand for the highest "experience" 
level in each billet. More precisely, using the "ascending ladder 
height," 

r(t) = (r>) r* (0) 

1 K* 

where 

r* (t) = r**(t) for u = 1, 2, . . ., K* 
u Su 

with r’^*(0) = r”^* as given in the previous section. 

Step 2: Update r*(t) by the usual Markovian equation of 

r’"’"(t+l) = r’"(t)P* 

steps: Update r**(t) by "shifting" its components one component to 

the right, that is, 

0 if j = 1 

"■j ifj = 2 

with r**(0) = r’^* as in Step 1. 
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step 4: Replace those components of the updated vector, r**(t+l), which 

stand for zero "experience" level in each billet by successive 
components of the updated vector r*(t+l). More precisely, 

r** (t+l) = r* (t+1) for u = 0, 

Su+1 u+1 



Step 5: As a last step in the updating process the accession vector, c** is 

added to the vector r_**(t+l) obtained in Step 4. That is, 

new r**(t+l) = old r**(t+l) + c ** 

These five steps are now repeated for the next higher t-value until r**(T) 
is obtained. Then Steps 6, 7, and 8 follow: 

Step 6: To add the appropriate components of r *’^(T) is easiest 

accomplished in terms of the "ascending ladder heights," Su- 
Namely, define 

a*(T) = (q^’'(T), ..., q^^(T)) for T = l,2,... 



where 

q* ,(T)= X for u = 0, 1, ..., K*-l 

^ + 1 j=Su+l j 

Step 7: The vector g*(T) is expanded to insert zeroes at the appropriate 

places, representing infeasible billets. This is accomplished in 
terms of the "ascending ladder epochs," kj. Define 

g(T) = (q^(T) q^(T)) 

where 
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when k = kj for j = 1, K* 



q,m = 



-< 



q (T) 

i 

0 



otherwise. 



Step 8: Finally, the q(T) is re-raveled into the original data format of an 

A by R matrix. That is, the matrix Q(T) is defined with elements 

where k = (n-l)A+i, with i = 1, ..., A and n = 1, ..., R. 

The numerical example introduced at the end of the previous section 
is continued to illustrate the procedure outlined above. The forecasting 
horizon of T=2 will be used. 

In Step 1, the "ascending ladder heights" 

So = 0, Si = 2, S2 = 5, S3 = 6, S4 = 8, S5 = 11, S6 = 13, S7 = 17 
are used to find the K* = 7 dimensional vector 

r(0) = (25, 30, 10, 15, 20, 25, 20). 

In Step 2, Markovian forecasting is applied to this vector using the 7 by 7 
P*-matrix given at the very end of the previous section. Thus, 
r*(l) = r*(0)P* = (0, 0, 36.5, 13, 18.5, 12, 2). 

In Step 3, the K** = 17 dimensional vector, r*’^(0) = r** is shifted one 
component to the right resulting in 

r**(l) = (0, 25, 25, 30, 30, 30, 10, 15, 15, 20, 20, 20, 25, 25, 20, 20, 20). 

Next, in Step 4, using the "ascending ladder heights," again, appropriate 
components of r**(l) above are replaced by successive components of r**(l) 
computed in Step 2. The result is 
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r**(l) = (0, 25, 0, 30, 30, 36.5, 13, 15, 18.5, 20, 20, 12, 25, 2, 20, 20, 20) 

where the replaced components are marked by an arrow above them. 
In Step 5, the accession vector 

c*”^ = (70, 0, 50, 0, 0, 0, 0, 0, 10, 0, 0, 0, 0, 0, 0, 0, 0) 



is added to the result of Step 4 to obtain 



r**(l) = (70, 25, 50, 30, 30, 36.5, 13, 15, 28.5, 20, 20, 12, 25, 2, 20, 20, 20). 



Next, Steps 1 through 5 are repated for t = 1 with the following results: 
Step 1 : r»(l) = (25, 30, 36.5, 1 5, 20, 25, 20) 

Step 2: r(2) = r*(l)P* = (0, 0, 36.5, 13, 39.7, 12, 2) 

Step 3: r»»(2) = (0, 70, 25, 50, 30, 30, 36.5, 13, 15, 28.5, 20, 20, 12, 25, 2, 20, 20) 

Step 4: r*(2) = (0, 70, 0, 50, 30, 36.5, 13, 13, 39.7, 28.5, 20, 12, 12, 2, 2, 20, 20) 



Step 5: new r **(2) = (70, 70, 50, 50, 30, 36.5, 13, 13, 49.7, 28.5, 20, 12, 12, 2, 2, 20, 20) 
Since r**(T) has been obtained, it remains only to reformat the result. 

In Step 6, the "ascending ladder heights," are used again to compute 
components of the vector 

3*(2) = (140, 130, 36.5, 26, 98.2, 24, 44). 

The details are as follows: 




q*(2)= Sr**(2)= 49.7 +28.5 +20 = 98.2 q*(2)= Z r**(2) = 12 + 12 = 24 

^ j=9 j ^ j=12 i 



20 



and 



17 

q!(2) = X = 2 + 2 + 20 + 20 = 44. 

‘ j=14 j 

In Step 7, the "ascending ladder epochs," 

ki = 1, lc2 = 3, k3 = 5, k4 = 6, ks = 7, k6 = 1 
are used to expand the K* = 7 dimensional vector 3.=^(2) into the K=12 
dimensional vector 

A 

^(2) = (140, 0, 130, 0, 36.5, 26, 98.2, 0, 0, 0, 24, 44). 

Finally, in Step 8, this vector is re-raveled into the 3 by 4 matrix 







140 


0 


98.2 


0 




Q(2) = 




0 


36.5 


0 


24 






\ 


130 


26 


0 


44 


y 



In practice, each element of this matrix would likely be rounded to the 
nearest integer if this is indeed the end of the forecasting process, as these are 
the expected numbers of officers two quarters into the future. 

4. AN APPLICATION. 

A typical application of the model to the Surface Warfare Officer (SWO) 
community of the U.S. Navy follows. The choice of activity types, listed in 
Table 1, was dictated by both the nature of the community and the topical 
issue of "jointness" referred to in the Introduction and Background. Thus, to 
analyze the influence of the new legislation on SWO career paths the activity 
types of "JPME" (Joint Professional Military Education) and "Joint Tours," 
were broken out separately in addition to "Postgraduate Education" and 
"SWO Education" (such as Department Head School) from the other "Shore 
Duty," with "Fleet Unit" being the main activity why the community exists. 
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The tourlength matrix, L, the Incumbents matrix, Q, and the Accessions 
matrix, C, are given in Tables 2, 3, and 4, respectively. Finally, the eleven 
transition probability matrices are presented in Table 5. 

The model was then used to forecast the distribution of SWOs four, ten 
and forty quarters into the future with Tables 6, 7 and 8 showing the results. 
The CPU-times required to perform these calculations were, respectively, 7.69, 
16.31, and 59.1 seconds on a Zenith Z-248(R) microcomputer with 1.1 MB 
RAM. 

The model's user-friendly features which enable the user to effect changes 
in any of the data input values given in Tables 1 through 5 and to carry out a 
billets versus officers comparison in a reasonably efficient manner are 
detailed in Johnson [3] and are not discussed here. 
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TABLE 1. CURRENT LIST OF ACTIVITY NAMES 
IN THE SURFACE COMMUNITY 



1. POSTGREDUC 

2. JPME 

3. JOINT TOUR 

4. SWOEDUC 

5. FLEET UNIT 

6. SHORE DUTY 



TABLE 2. CURRENT TOUR LENGTHS IN QUARTERS 



ACnVITIES/TOURS 


1 


2 


3 


4 


5 


6 


7 


8 


9 


ID 


u 


12 


1 . 


POSTGR EDUC 


0 


8 


8 


8 


6 


6 


6 


6 


6 


0 


6 


0 


2. 


JPME 


0 


0 


2 


0 


2 


2 


2 


2 


2 


2 


2 


2 


3. 


JOINT TOUR 


0 


8 


8 


8 


8 


8 


8 


8 


8 


8 


8 


8 


4. 


SWO EDUC 


0 


2 


2 


2 


2 


2 


0 


0 


0 


0 


0 


0 


5. 


FLEET UNIT 


12 


6 


6 


6 


6 


6 


6 


6 


6 


9 


9 


9 


6. 


SHORE DUTY 


0 


8 


8 


8 


8 


8 


8 


8 


8 


8 


8 


8 



TABLE 3. CURRENT INCUMBENTS 



AcnvixiEsyrouRS 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1. 


POSTGR EDUC 


0 


124 


133 


12 


12 


28 


8 


7 


2 


0 


2 


0 


2. 


JPME 


0 


0 


3 


0 


2 


11 


9 


12 


15 


7 


4 


1 


3. 


JOINT TOUR 


0 


1 


3 


3 


3 


12 


15 


36 


91 


41 


39 


19 


4. 


SWO EDUC 


0 


24 


55 


69 


15 


9 


0 


0 


0 


0 


0 


0 


5. 


FLEET UNIT 


3744 


765 


168 


299 


465 


253 


225 


194 


313 


172 


69 


46 


6. 


SHORE DUTY 


0 


466 


530 


126 


88 


119 


117 


244 


428 


252 


137 


83 



TABLE 4. CURRENT ACCESSIONS 



ACTIVITIESyTOURS 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


n 


12 


1 . 


POSTGR EDUC 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2. 


JPME 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


3. 


JOINT TOUR 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


4. 


SWO EDUC 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


5. 


FLEET UNIT 


315 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


6. 


SHORE DUTY 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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TABLE 5. CURRENT TRANSFER RATES 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 1 



AcnvrnEs/ACTrviTOs 


1 


2 


3 


4 


5 


6 


1. POSTGR EDUC 


0 


0 


0 


0 


0 


0 


2 JPME 


0 


0 


0 


0 


0 


0 


a JOINT TOUR 


0 


0 


0 


0 


0 


0 


4 SWOEDUC 


0 


0 


0 


0 


0 


0 


S FLEET UNIT 


ao6 


0 


0 


0B6 


040 


0J8 


shore duty 


0 


0 


0 


0 


0 


0 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 2 



AcnvrnEs/Acnvrrres 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


uoo 


0 


0 


2 JPME 


0 


0 


0 


0 


0 


0 


a JOINT TOUR 


0 


0 


0 


0 


0 


0 


4 SWO EDUC 


0 


0 


0 


0 


uoo 


0 


S FLEET UNIT 


0.13 


0 


0 


0 


OJO 


052 


6. SHORE DUTY 


0 


0 


0 


024 


0 


0 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 3 



ACT! viTi Es/ Acrrv rriEs 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


IBO 


0 


0 


2 JPME 


0 


0 


0 


0 


0 


0 


a JOINT TOUR 


0 


0 


0 


0 


0 


0 


4 SWO EDUC 


0 


0 


0 


0 


uoo 


0 


S FLEET UNIT 


0 


0 


0 


0 


uoo 


0 


k SHORE DUTY 


0 


0 


0 


053 


0 


025 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 4 



ACT! vm Es/Acnv rriEs 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


0 


0 


0 


2 JPME 


0 


0 


0 


0 


0 


0 


a JOINT TOUR 


0 


0 


0 


0 


0 


0 


4 SWO EDUC 


0 


0 


0 


0 


uoo 


0 


S FLEET UNIT 


005 


0 


0 


0 


050 


0.15 


6. SHORE DUTY 


0 


0 


0 


050 


0 


0 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 5 



ACT! VITIES/Acnv rriES 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


0 


lOO 


0 


2 JPME 


0 


0 


0 


0 


0 


0 


a JOINT TOUR 


0 


0 


0 


0 


0 


0 


4 SWO EDUC 


0 


0 


0 


0 


lOO 


0 


S FLEET UNIT 


005 


007 


002 


0 


056 


040 


k SHORE DUTY 


0 


0 


0 


0 


050 


0 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 7 



AcnvmEs/AcnvrriEs 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


0 


lOO 


0 


2 JPME 


0 


0 


050 


0 


050 


0 


a JCXNTTOUR 


0 


0 


0 


0 


lOO 


0 


4 SWO EDUC 


0 


0 


0 


0 


0 


0 


S FLEET UNIT 


006 


0.15 


006 


0 


025 


050 


k SHORE DUTY 


0 


0 


0 


0 


lOO 


0 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 8 



ACT! vrnEs/AcnvmEs 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


0 


lOO 


0 


2 JPME 


0 


0 


0.90 


0 


050 


0 


a JOINT TOUR 


0 


0 


0 


0 


lOO 


0 


4 SWO EDUC 


0 


0 


0 


0 


0 


0 


5 FLEET UNIT 


0 


0.15 


0.15 


0 


025 


045 


k SHORE DUTY 


0 


0 


0 


0 


Oj65 


055 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 9 



AcnvmEs/ACTivrrres 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


0 


0 


0 


2 JPME 


0 


0 


050 


0 


050 


0 


a JOINT TOUR 


0 


0 


0 


0 


lOO 


0 


4 SWO EDUC 


0 


0 


0 


0 


0 


0 


& FLEET UNIT 


0 


006 


002 


0 


0.15 


0j65 


k SHORE DUTY 


0 


0 


0 


0 


45 


025 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 10 



ACT! VmES/ACnV [TIES 


1 


2 


3 


4 


5 


6 


1. POSTGR EDUC 


0 


0 


0 


0 


0 


0 


2 JPME 


0 


0 


050 


0 


050 


0 


a JCftNT TOUR 


0 


0 


0 


0 


0.95 


0 


4 SWO EDUC 


0 


0 


0 


0 


0 


0 


5 FLEET UNIT 


0 


0 


0.10 


0 


0.15 


040 


k SHORE DUTY 


0 


0 


0 


0 


020 


050 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 11 



AcnvrnES/AcnviTiES 


1 


2 


3 


4 


5 


6 


1. POSTGR EDUC 


0 


0 


0 


0 


0 


0 


2 JPME 


0 


0 


0 


0 


0 


0 


a JOINT TOUR 


0 


0 


0 


0 


050 


0 


4 SWO EDUC 


0 


0 


0 


0 


0 


0 


5 FLEET UNIT 


0 


0 


020 


0 


0 


055 


6. SHORE DUTY 


0 


0 


0 


0 


050 


055 



CURRENT TRANSFER RATES WHEN 
LEAVING TOUR NUMBER 6 



A cn vrri ES/ACTIVITIES 


1 


2 


3 


4 


5 


6 


L POSTGR EDUC 


0 


0 


0 


0 


UOO 


0 


2 JPME 


0 


0 


040 


0 


0.60 


0 


a JOINT TOUR 


0 


0 


0 


0 


055 


0 


4 SWO EDUC 


0 


0 


0 


0 


0 


0 


S FLEET UNIT 


ox» 


OJO 


OJM 


0 


040 


040 


k SHORE DUTY 


0 


0 


0 


0 


055 


0 
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TABLE 6. EXPECTED NUMBERS OF OFFICERS 4 QUARTERS FROM NOW 



ACnVITIES/TOURS 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 . 


POSTGR EDUC 


0 


126 


135 


8 


14 


26 


10 


10 


0 


0 


0 


0 


2. 


JPME 


0 


0 


0 


0 


0 


11 


8 


11 


10 


6 


0 


0 


3. 


JOINT TOUR 


0 


0 


0 


0 


0 


14 


24 


37 


84 


37 


35 


14 


4. 


SWO EDUC 


0 


31 


60 


78 


16 


0 


0 


0 


0 


0 


0 


0 


5. 


FLEET UNIT 


3756 


755 


162 


328 


464 


272 


235 


183 


243 


280 


103 


61 


6. 


SHORE DUTY 


0 


457 


530 


130 


74 


185 


127 


200 


317 


317 


137 


75 



TABLE 7. EXPECTED NUMBERS OF OFFICERS 10 QUARTERS FROM NOW 



ACnVITIES/TOURS 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1. 


POSTGR EDUC 


0 


125 


131 


0 


16 


23 


14 


12 


0 


0 


0 


0 


2. 


JPME 


0 


0 


0 


0 


0 


11 


9 


12 


9 


4 


0 


0 


3. 


JOINT TOUR 


0 


0 


0 


0 


0 


12 


32 


33 


79 


26 


29 


14 


4. 


SWO EDUC 


0 


31 


58 


77 


17 


0 


0 


0 


0 


0 


0 


0 


5. 


FLEET UNIT 


3774 


749 


169 


342 


495 


276 


267 


192 


193 


373 


150 


80 


6. 


SHORE DUTY 


0 


449 


526 


132 


64 


248 


142 


155 


190 


319 


158 


73 



TABLE 8. EXPECTED NUMBERS OF OFFICERS 40 QUARTERS FROM NOW 



ACnVITIES/TOURS 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1. 


POSTGR EDUC 


0 


126 


131 


0 


17 


25 


14 


17 


0 


0 


0 


0 


2. 


JPME 


0 


0 


0 


0 


0 


12 


10 


17 


13 


4 


0 


0 


3. 


JOINT TOUR 


0 


0 


0 


0 


0 


13 


34 


42 


113 


31 


34 


20 


4. 


SWO EDUC 


0 


32 


59 


76 


16 


0 


0 


0 


0 


0 


0 


0 


5. 


FLEET UNIT 


3780 


756 


170 


346 


505 


289 


340 


255 


229 


311 


132 


73 


7. 


SHORE DUTY 


0 


454 


524 


131 


69 


269 


153 


227 


232 


243 


166 


84 
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